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Abstract — Chaotic sequences have been widely used as 
pseudorandom sequences. But up to now, we have no simple and 
effective method to test and improve their autocorrelation and RFM 
autocorrelation performances. In the paper, by method of phase 
space, relations between autocorrelation and phase space trajectory 
were studied, simple and effective methods were presented to assess 
and improve their autocorrelation and RFM autocorrelation 
performances. Many simulations were presented, too. 
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I. INTRODUCTION 

Chaotic sequence has been widely used in communication and 
radar as pseudorandom sequence [125] . But we have no simple 
method to test and improve their autocorrelation and modulated 
autocorrelation performances. Therefore, applications of chaotic 
sequences are limited [118] . (For easy to describe, in the paper, 
autocorrelation of modulated signal is referred to as modulated 
autocorrelation.) 

Recently, chaotic sequences were used in noise radar as noise 
source [111,18] . But because of structure of chaotic attractor, some 
chaotic sequences have poor autocorrelation performance, or 
after modulation, such as Random Frequency Modulation 
(RFM) [36] , autocorrelation performance of modulated signal is 
poor [36] . This problem impede the detection and identification of 
chaos-based signals, so autocorrelation and modulated 
autocorrelation of chaotic sequence have attracted many 
academic attention [27] . 
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In the paper, by phase space, we proved some theorems to 
present simple and effective methods for assessment and 
improvement of autocorrelation and RFM autocorrelation 
performances of chaotic sequences. Simulations were presented 
to verify them. Too. 



II. The problem of chaotic sequence autocorrelation 
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(a) Normalized autocorrelation 
of Tent sequence 



(b)Normalized autocorrelation of 
RFM modulated Tent sequence 



Fig. 1. Normalized autocorrelation and RFM modulated autocorrelation 
functions of Tent sequence 



-5 

n -10 

^-15 

^-20 

-25 

-3Q 




■5-4-3-2-1012345 

m 
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of Bernoulli sequence RFM modulated Bernoulli sequence 

Fig. 2. Normalized autocorrelation and RFM modulated autocorrelation 
functions of Bernoulli sequence 

For communication, the shape of a good autocorrelation 
function ought to like a sharp needle, as shown in figure 1(a). It 
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has a sharp and outstanding main peak, and has no outstanding 
side lobe, so this signal can be easily detected and identified. 

The autocorrelation function of Tent sequence is presented in 
figure 1(a), its main peak is sharp and outstanding, and its side 
lobe is low and can hardly be find out, so its performance is very 
good. Although Tent sequence autocorrelation is good, but after 
Random Frequency Modulation (RFM) [36] , its RFM 
autocorrelation performance becomes bad, as shown in figure 
1(b), there are many outstanding side lobes. 

The autocorrelation function of Bernoulli sequence is shown in 
figure 2(a), its main peak is blunt and fat, so its performance is 
poor. But after RFM, its RFM autocorrelation performance 
becomes good, as shown in figure 2(b). 

Therefore, some chaotic sequences have bad autocorrelation 
performance, and some chaotic sequences have poor modulated 
autocorrelation performance. So it needs to study this 
phenomenon, and find out its rule to assess and improve their 
autocorrelation and RFM autocorrelation performances for better 
applications. 



III. Autocorrelation of chaotic sequence 

By Autocorrelation Phase-space Axis Symmetric(APAS) 
theorem [25] , we can estimate autocorrelation performance of a 
sequence. It only needs to test whether its trajectory in delay 1 
phase space is x-axis or y-axis symmetrical. If its trajectory is 
x-axis or y-axis symmetrical, then its autocorrelation 
performance is good, namely, its autocorrelation has sharp and 
outstanding main peak, and has no remarkable sidelobe. 

we list APAS theorem as follows. 
Theorem 1: A stationary ergodic discrete real dynamic system 
x (n + 1) = <p[x(n)] (1) 

Where (p(x) is a single valued function, value area of x(n) is [-a, 

a], and 'a' is a positive real number; the mean of {x(n)} is zero, 
values of x(n) are statistically balanced, namely 
f[x(n)] = f[-x(n)] , where, f[x(n)] is probability density 
function of x(n). 

If (p(x) satisfy any one of following condition 

x{n + 1) = <p[x(n)] = (p[-x(n)] (2) 

f[x(n\ x(n + 1)] = f[-x(n), x(n + 1)] (3) 
Then, we get following conclusion 

when N — > °o , and delay m ^ , to any other delay m , the 
normalized autocorrelation function of {x(n)} is near to zero. 
That is 



lim R{m) = lim V x(n)x(n + m) = 

*->- ^-r(0)ti ' (4) 

me N, m^O 

Theorem 2: A stationary ergodic discrete real dynamic system 
x(n + 1) = <p[x(n)] (5) 

Where inverse function (p~ (x) of (p(x) is single valued. 

Value area of x(n) is [-a, a], and 'a' is positive real; the mean of 
{x(n)} is zero, values of x(n) are statistically balanced, namely 
f[x(n)] = f[-x(n)] , where, f[x(n)] is probability density 
function of x(n). 

If (p(x) satisfy any one of following condition 

x(n) = (p~ l [x(n + 1)] = q>~ 1 [-x(n + 1)] 
f[x(ri), x(n + 1)] = f[x(n) -x(n + 1)] 
Then, we get following conclusion 

when N — > oo 9 and delay m ^ , to any other delay m , the 
normalized autocorrelation function of {x(n)} is near to zero. 
That is 

1 A/-1 

lim R(m) = lim V x(n)x(n + m) = 

A^oo A^oo r (()) ^ 

meN, m*0 



(6) 
(7) 



Combining Theorem 1 and Theorem 2, we call them as 
Autocorrelation Phase-space Axis Symmetric(APAS) Theorem. 

APAS theorem was proved and verified in paper[25]. 

From APAS theorem we can see, to estimate a sequence 
autocorrelation performance, it only needs to test whether its 
delay 1 phase space trajectory is x-axis or y-axis symmetrical. If 
its trajectory is x-axis or y-axis symmetrical, then its 
autocorrelation performance is good, and has sharp and 
outstanding main peak, and has no remarkable sidelobe. 

We have assessed and improved autocorrelation performances 
of some chaotic sequences [26 ' 27] by APAS theorem, it is verified 
that APAS theorem is simple and effective. 

IV. RFM Autocorrelation of chaotic sequence 

Recently, Random Frequency Modulation(RFM) received 
more and more academic attention [37] . 

A discrete chaotic system, its dynamic equation is 

x(n + 1) = <p[x(n)] (9) 

Value area of x(n) is [-a, a], and 'a' is maximum of x(n), is a 
positive real number; the mean of {x(n)} is zero, values of x(n) 
Sire statistically balanced, namely f[x(n)] = f[— x(n)] , where, 

f\x(n)\ 

JL v n is probability density function of x(n). 

Discrete RFM signal of {x(n)} with complex envelope [3 ' 4] 
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s(i) = A exp[ jlxKJ^ x(n) + j27T<p ] 



(10) 



Where K is an index which combine modulation index and 
sample rate, is a positive real number. 

For satisfy the Nyquist criterion, must meet condition of 

2Ka <1 (11) 

We choose 

2Ka = 1 (12) 

Then, phase angle of ^./ ^ ^ is distributed in a 

single full cycle [-71, k] . In order to convenient for description, we 
call equation (12) as Single Full Cycle(SFC) condition. Without 
loss of generality, we choose k=l, a=0.5, and following 
simulations all satisfy the condition. 

A stationary ergodic sequence {x(n)| with length N, its RFM 
autocorrelation function r(m) is 



N-m 



r(m)= 2^s(n)s*(n + m) 



n=\ 



Normalized RFM autocorrelation function R(m) is 

i N-m 

Rirn) = — - ^s(n)s*(n + m) 



(13) 
(14) 



RFM-APTS Theorem: A stationary ergodic discrete real 
dynamic system 

x(n + 1) = <p[x(n)] (15) 

Value area of x(n) is [-a, a], the mean of {x(n)} is zero, values of 
x(n) are statistically balanced, namely f[x(n)] = f[— x(n)] , 

where, f[x(n)] is probability density function of x(n). RFM 

signal of {x(n)| satisfy SFC condition. 

On delay 1 trajectory of {x(n)}, if the trajectory is statistically 
translational symmetrical, as shown in figure 3, namely 



x(n + 1) = q>[x(n)\ = q>[x(n) + a] 

x(n) < 

x(n + 1) = <p[x(n)] = (p[x(n) - a] 

x(n) > 



(16) 



and 



f[x(n)] = f[x(n) + a] 



f[x(n)] = f[x(n)-a] 



x(n) < 



x(n) > 



(17) 



Where, /[x(n)] is probability density function of x(n) 

Then, we get following conclusion 

when N — > °° , and delay m ^ , to any other delay m , the 
normalized RFM autocorrelation function of {x(n)} is near to 
zero. That is 

i AM 

lim R(m) = lim V s(n)s* (n + m) « 

AM- ^ ^ N ^ r ^^ ^ ^ ^ (1g) 

me N, m^O 
Proof: 
(1) prove lim R(l) = . 

AM- 

On delay 1 phase trajectory, we can arbitrary choose a couple 
of points D(k) , D(j) which are translational symmetrical. D(k) is 
on right part of phase plane, and coordinate is [xi(k), Xi(k+l)] ; D(j) 
is on left part of phase plane, and coordinate is [x\(j), JCi(/+l)]. 
Then, have Xi(j)+a= x x (k) , Xi(k+l)= Xi(/+1), as shown in figure 
3. Figure 3 is an example figure using Bernoulli sequence as an 
illustration. 




Fig. 3. delay 1 phase trajectory of Bernoulli sequence 



Because 



f[x(n)] = f[x(n) + a] 



f[x(n)] = f[x(n)-a] 



x(n) < 



x(n) > 

Then, without loss of generality, choose x(n) > , 
so f[x(n)] = f[x(n)- a] o 

Let P[x(n + l) = x 1 (k + 1)| x(ri) = x x (k)] 

f[x(n + 1) = x x (k + 1)| x(n) = x x (k)] are conditional 
probability and conditional probability density function of 
x{n + Y) = x l (k+ 1) under condition x{n) - x 1 (k) 
respectively. In the paper, they are abbreviated as 
P[x { (k + 1)\ *,(*)], /[^(fc + l)! Xl (k)]. 
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To a point D(k)[xi(k), Xi(k+l)] on trajectory, because 
X x (k + 1) = (p[x l (k)] , yields 

PIXOfc + 1)! x l (k)] = P[x l (k + l)\ x l (k)-a] = l 

Then 

f[x x (k + \)\ x 1 (k)] = f[x 1 (k + l)\ x { (k)-a] 

Let f[x(n) = x l (&), x(n + l) = x l (k + 1)] is Joint 

probability distribution density function of x(n) = X x {k) and 



x(n + 1) = X x (k + 1) , and is abbreviated 



as 



f[ Xl (k), x,(fc + l)]. 

So, on trajectory, to translational symmetrical couple of points 
D(k), D(j), have 

f[x x (k\ x,(k + l)] 

=f[ Xl (k+i)\ x.mx/ix.m 

= f[x l (k + l)\ jc 1 (^)-a]x/[jc 1 (^)-a] 
= /[x 1 (j), ^O' + l)] 



(19) 



That is, distribution probability densities of the two points 
are identical, namely /[£>(&)]= /[£>(/)]. 

To a point D(k) [x^k), JCi(Jfc+l)], 
Let function h (k) is 

h(k) = s(k-l)s\k) 

= A 2 exp[-j2xKx(k)] (20) 

The couple of points D(k) > D(j), satisfy SFC condition, we 
get 



h(j) = A 2 exp[-j27iKx(j)] 
= A 2 exp[-j(27lKx(k) + 2xKa)] 
= -A 2 exp[-j27iKx(k)] 

= -h(k) 
Then, we get h(k) + h(j) = . 

On trajectory, to any near translational symmetrical two 
points D(k), D(j), have h(k) + h(j) « . 

Recall that translational symmetrical two points D(k)^ 
D(j) satisfy j[D(k)]= flPffl 



(21) 



Then, all the (N-l) points on the phase figure can be 
approximately divided into (N-l)/2 groups, and two points 
of each group are near translational symmetrical, and to 

each group D(k n ), D(j n ), satisfy h(k n ) + h(j n )~0. 

There are few points can't be grouped as they are not 
translational symmetrical or number of points is odd number. But 
distribution probability densities of any two points on trajectory 

are identical, when ™ ~~* °° , the distribution of points on 
trajectory tends to be more close to translational symmetrical, 
whereas number of points which can't be grouped is limited, and 
because of 

1 AM 

lim7?(l) = lim Ys(n)s\n+1) 

a 2 N-l 1 N-l 

= V exp \rj2jKin + 1)] = lim V Hn + 1) 

N-l 

Denominator lim r(0) = lim V I s(n) I 2 — > +°o 

n=l 

So, the effects of no grouped points to R(l) tend to be zero, 
We get 

i N-l 

lim/?(l)=lim Yh(n + 1) 



1 



(N-l)/2 



lim Y0~0 



1 n=l 

(2) prove lim R(m) « 

N^oo 

Based on lim R(l) ~ , Similar to the proof of APAS 

N^oo 

Theorem [25] , we can obtain lim R{m) « . 
Which completes the proof. 

V. Conclusions about assessment of Autocorrelation 
and RFM Autocorrelation 

Sum up, we get conclusions that, 

(l)Assessment of autocorrelation: 

A stationary ergodic statistically balanced time sequence 
which mean is zero, if its delay 1 phase trajectory is statistically 
axial symmetrical, 

Then, its autocorrelation performance are good. 

(2)Assessment of RFM autocorrelation: 

A stationary ergodic statistically balanced time sequence 
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which mean is zero and satisfy SFC condition, if its delay 1 phase 
trajectory is statistically translational symmetrical, 

Then, its RFM performance are good. 

From APAS and RFM-APTS theorems we can see, to estimate 
a sequence autocorrelation and RFM autocorrelation 
performances, only needs to test the shape of delay 1 phase space 
trajectory, if it is x-axis or y-axis symmetrical, then 
autocorrelation performance is good; if it is translational 
symmetrical, then RFM autocorrelation performance is good. 

Phase space trajectory figure is an effective tool to observe a 
dynamic system internal characteristic, by which we can 
recognize structure of the system. Generally, compared with other 
larger delay phase space figure, the structure of delay 1 phase 
space figure is far more simple and clear, so APAS and 
RFM-APTS theorems is a simple and effective method to judge a 
sequence autocorrelation and RFM autocorrelation 
performances. 

VI. Testing autocorrelation and RFM autocorrelation 

PERFORMANCES BY APAS AND RFM-APTS THEOREMS 

Tent sequence phase trajectory is shown in figure 4(a), Its 
trajectory is y-axis symmetrical, satisfies axial symmetrical 
condition, but doesn't satisfy translational symmetrical condition. 
By APAS and RFM-APTS theorems, its autocorrelation 
performance is good, RFM autocorrelation performance is poor. 
Its normalized autocorrelation and RFM autocorrelations are 
shown in figure 4(b), (c), and meet APAS and RFM-APTS 
theorems. 

Delay 1 phase trajectory of Bernoulli sequence is shown in 
figure 5(a), its trajectory doesn't satisfy axial symmetrical 
condition, but satisfies translational symmetrical condition. By 
APAS and RFM-APTS theorems, its autocorrelation 
performance is poor, RFM autocorrelation performance is good. 
Its normalized autocorrelation and RFM autocorrelations are 
shown in figure 5(b), (c), and meet APAS and RFM-APTS 
theorems. 




x(n) 
(a) Delay 1 phase space trajectory of Tent sequence 
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(a) Normalized autocorrelation 
of Tent sequence 



(b)Normalized RFM autocorrelation 
of Tent sequence 



Fig.4. Delay 1 phase space trajectory of Tent sequence and 
its autocorrelation and RAM, RFM, RPM autocorrelation functions 
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(a) Delay 1 phase space trajectory of 
Bernoulli sequence 
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(b) Normalized autocorrelation 
of Bernoulli sequence 



(c)Normalized RFM autocorrelation 
of Bernoulli sequence 

Fig. 5. Delay 1 phase space trajectory of Bernoulli sequence and 
its autocorrelation and RAM, RFM, RPM autocorrelation functions 
Skew Tent sequence is different from upper Tent and Bernoulli 
sequences. Delay 1 phase trajectory of Skew Tent sequence is 
shown in figure 6(a), its trajectory neither satisfies axial 
symmetrical condition, nor satisfies translational symmetrical 
condition. By APAS and RFM-APTS theorems, its 
autocorrelation and RFM autocorrelation performances all are 
poor. Its normalized autocorrelation and RFM autocorrelations 
are shown in figure 6(b), (c), and meet APAS and RFM-APTS 
theorems. 
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(a) Delay 1 phase space trajectory 
of Skew Tent map which vertex is 0.4 
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(b) Normalized autocorrelation 
of Skew Tent sequence 
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(c)Normalized RFM autocorrelation 
of Skew Tent sequence 



Fig. 6. Delay 1 phase space trajectory of Skew Tent sequence and 
its autocorrelation and RAM, RFM, RPM autocorrelation functions 

VII. Improvement of autocorrelation and RFM 

AUTOCORRELATION PERFORMANCES BY APAS AND RFM-APTS 
THEOREMS 

To a chaotic map which trajectory is not axial or translational 
symmetrical, that doesn't satisfy conditions of APAS and 
RFM-APTS theorems, we can modify its map, by methods of 
compressing, copy, fold, mirror, translation, etc, change its 
attractor to satisfy axial and translational symmetrical conditions, 
then the modified chaotic sequences have good autocorrelation 
and RFM autocorrelation performances. 

To tent map, its trajectory is axial symmetrical, so we modify its 
map, copy, compress and translate its trajectory, we call the 
modified map as Axial and Translational Symmetrical 
Tent(ATS-Tent) map. Phase trajectory of ATS-Tent is shown in 
figure 7(a), is Axial and Translational Symmetrical, therefore its 
autocorrelation and RFM autocorrelation performances are good, 
as shown in figure 7(b),(c), and meet APAS and RFM-APTS 
theorems. 

Or 

-5 
ffl -10 
^-15 
^-20 

-25 

-30, 



0.25 




-0.25 



(a) phase space trajectory of 
ATS-Tent sequence 
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(b) Normalized autocorrelation 
of ATS-Tent sequence 



■5-4-3-2-1012345 



(c) Normalized RFM autocorrelations of ATS-Tent sequence 



Fig.7. Delay 1 phase space trajectory of ATS-Tent sequence and 
its autocorrelation and RFM autocorrelation functions 

To Bernoulli map, we copy or mirror its trajectory two times, 
We call the map as Axial and Translational Symmetrical 
Bernoulli(ATS-Bernoulli) map, its delay 1 phase figure is shown 
in figure 8(a). It is Axial and Translational Symmetrical, then its 
autocorrelation and RFM autocorrelation performances all are 
good, as shown in figure 8(b), (c), and meet APAS and 
RFM-APTS theorems. 
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(a) phase space trajectory of 
ATS-Bernoulli sequence 
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(b) Normalized autocorrelation 
of ATS-Bernoulli sequence 




(c) Normalized RFM autocorrelation of ATS-Bernoulli sequence 

Fig. 8. Delay 1 phase space trajectory of ATS-Bernoulli sequence and 

its autocorrelation and RFM autocorrelation functions 

To Skew tent map, its trajectory neither is axial symmetrical, 

nor translational symmetrical. We copy or mirror its trajectory 

two times, We call the map as Axial and Translational 

Symmetrical Skew tent(ATS-Skew tent) map, its delay 1 phase 

figure is shown in figure 9(a). It is Axial and Translational 

Symmetrical, then its autocorrelation and RFM autocorrelation 

performances all are good, as shown in figure 9(b), (c), and meet 

APAS and RFM-APTS theorems. 





(a) phase space trajectory of 
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(b) Normalized autocorrelation of 
ATS- Skew Tent sequence 



CISME Vol.1 No. 2 2011 PP. 25-31 www.jcisme.org © World Academic Publishing 

30 



Communications in Information Science and Management Engineering 



CISME 



0, — ^ — ■ — - 
-5 
ffl -10 

■73 

^-15 

-25 

-30 ' — — ■ — - 
-5-4-3-2-1012345 

m 

(c) Normalized RFM autocorrelation of ATS- Skew Tent sequence 

Fig. 9. Delay 1 phase space trajectory of ATS- Skew Tent sequence and 

its autocorrelation and RFM autocorrelation functions 



VIII. Conclusion 

By phase space method, we presented APAS and RFM-APTS 
theorems for assessment and improvement of autocorrelation and 
RFM autocorrelation performances of chaotic sequences. It is a 
simple and effective method by APAS and RFM-APTS theorems 
to assess and improve performances of autocorrelation and RFM 
autocorrelation. Simulations were presented to verify the theorem 
and the method, too. 
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